Some methods of training radial basis neural networks in solving the Navier‐Stokes equations

The purpose of this research is to analyze the application of neural networks and specific features of training radial basis functions for solving 2‐dimensional Navier‐Stokes equations. The authors developed an algorithm for solving hydrodynamic equations with representation of their solution by the method of weighted residuals upon the general neural network approximation throughout the entire computational domain. The article deals with testing of the developed algorithm through solving the 2‐dimensional Navier‐Stokes equations. Artificial neural networks are widely used for solving problems of mathematical physics; however, their use for modeling of hydrodynamic problems is very limited. At the same time, the problem of hydrodynamic modeling can be solved through neural network modeling, and our study demonstrates an example of its solution. The choice of neural networks based on radial basis functions is due to the ease of implementation and organization of the training process, the accuracy of the approximations, and smoothness of solutions. Radial basis neural networks in the solution of differential equations in partial derivatives allow obtaining a sufficiently accurate solution with a relatively small size of the neural network model. The authors propose to consider the neural network as an approximation of the unknown solution of the equation. The Gaussian distribution is used as the activation function.     

On weak–strong uniqueness for the compressible Navier–Stokes system with non-monotone pressure law

We show the weak–strong uniqueness property for the compressible Navier–Stokes system with general non-monotone pressure law. A weak solution coincides with the strong solution emanating from the same initial data as long as the latter solution exists.     

On decay rate of solutions for the stationary Navier–Stokes equation in an exterior domain

In this paper, we consider the asymptotic behavior of an incompressible fluid around a bounded obstacle. By adapting the Schauder's estimate for stationary Navier–Stokes equation to improve the regularity, the problem is solved by using appropriate Carleman estimates. It should be noted that the minimal decaying rate for a general scalar equation is $\exp ?(?C|x{|}^{2+})$. However, the structure of the Navier–Stokes is special. Under the assumption for any nontrivial solution to be uniform bounded which is weaker than those in [10], we got the minimal decaying rate is $\exp ?(?C|x{|}^{\frac{3}{2}+})$ which is better than the results in general scalar cases.     

Some preliminary observations on a defect Navier–Stokes system

Some implications of the simplest accounting of defects of compatibility in the velocity field on the structure of the classical Navier–Stokes equations are explored, leading to connections between classical elasticity, the elastic theory of defects, plasticity theory, and classical fluid mechanics.     

Shape morphing and topology optimization of fluid channels by explicit boundary tracking

An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.     

The improvement of numerical modeling in the solution of incompressible viscous flow problems using finite element method based on spherical Hankel shape functions

In this paper, the finite element method with new spherical Hankel shape functions is developed for simulating 2‐dimensional incompressible viscous fluid problems. In order to approximate the hydrodynamic variables, the finite element method based on new shape functions is reformulated. The governing equations are the Navier‐Stokes equations solved by the finite element method with the classic Lagrange and spherical Hankel shape functions. The new shape functions are derived using the first and second kinds of Bessel functions. In addition, these functions have properties such as piecewise continuity. For the enrichment of Hankel radial basis functions, polynomial terms are added to the functional expansion that only employs spherical Hankel radial basis functions in the approximation. In addition, the participation of spherical Bessel function fields has enhanced the robustness and efficiency of the interpolation. To demonstrate the efficiency and accuracy of these shape functions, 4 benchmark tests in fluid mechanics are considered. Then, the present model results are compared with the classic finite element results and available analytical and numerical solutions. The results show that the proposed method, even with less number of elements, is more accurate than the classic finite element method.     

Accurate iteration-free mixed-stabilised formulation for laminar incompressible Navier–Stokes: Applications to fluid–structure interaction

Stabilised mixed velocity–pressure formulations are one of the widely-used finite element schemes for computing the numerical solutions of laminar incompressible Navier–Stokes. In these formulations, the Newton–Raphson scheme is employed to solve the nonlinearity in the convection term. One fundamental issue with this approach is the computational cost incurred in the Newton–Raphson iterations at every load/time step. In this paper, we present an iteration-free mixed finite element formulation for incompressible Navier–Stokes that preserves second-order temporal accuracy of the generalised-alpha and related schemes for both velocity and pressure fields. First, we demonstrate the second-order temporal accuracy using numerical convergence studies for an example with a manufactured solution. Later, we assess the accuracy and the computational benefits of the proposed scheme by studying the benchmark example of flow past a fixed circular cylinder. Towards showcasing the applicability of the proposed technique in a wider context, the inf–sup stable P2–P1 pair for the formulation without stabilisation is also considered. Finally, the resulting benefits of using the proposed scheme for fluid–structure interaction problems are illustrated using two benchmark examples in fluid-flexible structure interaction.     

Application of Reynolds stress transport turbulence closure models to flows affected by Lorentz and buoyancy forces

The effect of magnetic field strength and orientation on two types of electromagnetically influenced turbulent flows was studied numerically under the Reynolds averaged Navier–Stokes (RANS) framework. Previous work (Wilson et al., 2014) used an electromagnetically extended linear eddy-viscosity model, whilst the current paper focuses on the performance of a more advanced Reynolds stress transport type model both with and without electromagnetic modifications proposed by Kenjereš et al. (2004). First, a fully-developed 2D channel flow is considered with a magnetic field imposed in either the wall-normal or streamwise direction. Both forms of the RSM gave good agreement with the DNS data for the wall-normal magnetic field across the range of Hartmann numbers with the additional electromagnetic terms providing a small, but noticeable, difference. For the streamwise magnetic field, where electromagnetic influence is only through the turbulence, the electromagnetically extended RSM performed well at moderate Hartmann numbers but returned laminar flow at the highest Hartmann number considered, contrary to the DNS. The RSM results were, however, significantly better than the previous eddy-viscosity model predictions. The second case is that of unsteady 3D Rayleigh–Bénard convection with a magnetic field imposed in either a horizontal or vertical direction. Results revealed that a significant reorganization of the flow structures is predicted to occur. For a vertically oriented magnetic field, the plume structures increase in number and become thinner and elongated along the magnetic field lines, leading to an increase in thermal mixing within the core in agreement with Hanjalić and Kenjereš (2000). With a horizontal magnetic field, the structures become two-dimensional and a striking realignment of the roll cells’ axes with the magnetic field lines occurs. The results demonstrate the capability of the Reynolds stress transport approach in modelling MHD flows that are relevant to industry and offer potential for those wishing to control levels of turbulence, heat transfer or concentration without recourse to mechanical means.     

The ASBM-SA turbulence closure: Taking advantage of structure-based modeling in current engineering CFD codes

Structure-based turbulence models (SBM) carry information about the turbulence structure that is needed for the prediction of complex non-equilibrium flows. SBM have been successfully used to predict a number of canonical flows, yet their adoption rate in engineering practice has been relatively low, mainly because of their departure from standard closure formulations, which hinders easy implementation in existing codes. Here, we demonstrate the coupling between the Algebraic Structure-Based Model (ASBM) and the one-equation Spalart–Allmaras (SA) model, which provides an easy route to bringing structure information in engineering turbulence closures. As the ASBM requires correct predictions of two turbulence scales, which are not taken into account in the SA model, Bradshaw relations and numerical optimizations are used to provide the turbulent kinetic energy and dissipation rate. Attention is paid to the robustness and accuracy of the hybrid model, showing encouraging results for a number of simple test cases. An ASBM module in Fortran-90 is provided along with the present paper in order to facilitate the testing of the model by interested readers.